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Dynamics of two particles with short range repulsive or 
attractive interaction is studied numerically in the Harper 
model. It is shown that interaction leads to appearance of lo- 
calized states and pure-point spectrum component in the case 
when noninteractive system is quasi-diffusive or ballistic. In 
the localized phase interaction gives only stronger localization 
contrary to the case of two interacting particles in a random 
potential. 

PACS numbers: 05.45. +b, 72.15.Qm, 72.10.Bg 



The Harper model of electrons on a two-dimensional 
square lattice in the presence of a perpendicular magnetic 
field was intensively studied during last decades (see e.g 
). After fixing the quasimomentum in one of direc- 
tions the eigenvalue equation is reduced to a very simple 
form of one-dimensional quasiperiodic discrete chain 



2A cos(hn + (3)<f) n + <f>„ 



E(f> n 



(1) 



where the effective Plank's constant h/2ir gives the ratio 
of magnetic flux through the lattice cell to one flux quan- 
tum, j3 is a constant related to quasimomentum. For the 
original problem of electrons in a magnetic field the pa- 
rameter A should be fixed at A = 1 but generally one can 
consider the model (1) at different values of A. Intensive 
analytical and numerical studies [§|J^^] for typical irra- 
tional values Ti/2tt showed that for A > 1 the spectrum 
is a pure-point like with gapes and all eigenstates are ex- 
ponentially localized. For A < 1 the spectrum becomes 
continuous with extended eigenstates corresponding to 
ballistic classical motion. For A = 1 the situation is crit- 
ical with singular-continuous multifractal spectrum and 
power law localized eigenstates. 

While the properties of one-particle Harper model are 
now well understood in many respects the question about 
effects of particles interaction was not much investigated 
up to now. Indeed, the main physical problem of interact- 
ing particles at finite particle density is very complicated 
both for analytical and numerical investigations. One of 
the approaches to understanding of this problem is to 
analyze the effect of interaction of only two particles in 
the Harper model. Recently such approach has given a 
number of interesting results for two interacting particles 
(TIP) in a random potential showing that even repulsive 
particles can form a pair which can propagate on a large 
distance [^|J8) . In this Letter I address the problem of TIP 
in a quasiperiodic potential showing that here interaction 
effects can be quite different from the case of random 



potential. The investigation of such type of model will 
allow to analyze the stability of multifractal spectrum in 
respect to interactions. 

For two particles on a square lattice [x, y) with mag- 
netic flux and on site interparticle interaction the eigen- 
value equation has the form: 



Yxi+l,yi,X2,V2 
1 Px 1 ,yi + l,x 2 ,y2 



2/1 ,2:2 + 1,2/2 



T c YX\ — 1,1/1, X2, V2' 

+ V^i, i/i — l, X2,y 2 J t~ 
T c YX\,yx,X2 — 1,V2~ 



(2) 



' > Pxi,yi,X2, 3/2+1 + rxi,Ki,X2,3/2- 
J X! ,x 2 fiyi ,2/2 V^l ,2/1 ,12 ,J/2 



Udxi .xn&in At2 V^i ,2/1 ,^2 ,2/2 — E / ^Px±,yi,X2,y2 



Here (x, y) are integers marking the sites of the square 
lattice, the indices 1,2 note two particles, U is on site 
interaction, and Ti — 2w^)/<f>Q is determined by the ratio 
of magnetic flux <f> through the unit cell to the quantum 
of flux cf>Q. The direct investigation of equations (2) is a 
quite complicated problem. Therefore, I reduce it to a 
simpler one with Bloch waves propagating in x— direction 

^x l , yi ,x2.. V 2 = (p Vl , y2 J dkidk 2 A kuk2 exp(i(k 1 xi + k 2 x 2 )) 
that leads to TIP in the Harper model with effectively 
renormalized interparticle interaction: 

(2Acos(?ini + (3i) + 2\cos(hn 2 + (3 2 ) + US nun2 )ip nun2 + 

9-^1 + 1,712 ~t~ ^Pn\—l,n2 9^1,1,712 + 1 *Pni,Tl2 — 1 — ^^Pfll, fl2 

(3) 

Here for generality there is introduced the parameter A 
which should be taken equal to 1 for the model (2) and 
the coordinates y\^ 2 are replaced by n\^ 2 . The parame- 
ters Pi j2 are /?i,2 = fci,2 and the strength of renormalized 
interaction is U = U J dkAk 1 +k 2 -k,k/Ak 1 ,k 2 - The inves- 
tigation of the properties of the model (3) should allow 
to understand better the properties of the original TIP 
problem on 2d-lattice (2). 

The time evolution of models (2), (3) is governed by 
eqs. (2), (3) with E in the right hand side replaced 
by id/dt. This evolution for the model (3) was stud- 
ied numerically on effective 2-dimensional lattice with 
size N x N < 301 x 301. The flux ratio was fixed 
at the golden mean value h/2ir = (^5 - l)/2. For 
A = 1 the spectrum of noninteracting problem is singular- 
continuous with multifractal properties (^-|5). Due to 
that the spreading over the lattice is similar to the dif- 
fusive case with the second moments of probability dis- 
tribution a± =< («i ± n 2 ) 2 > growing approximately 
linearly with time (see Fig.l). The switched on inter- 
action leads to a significant decrease in the rate of this 
growth, namely approximately in 10 times for the case 
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in Fig.l. Here, initially at time t = two particles are 
located at the same site so that all probability is con- 
centrated at n\ = ri2 = (the same initial conditions 
were used in Figs. 2-3). The analysis of the probabil- 
ity distribution P± = J2 n ±=const p ( n ii n 2) dependence 
on n± — \n\ ± tt^I/v^ shows that its tail has Gaus- 
sian shape. However, while in the noninteractive case all 
the probabilities spread diffusively over the lattice in the 
case with interaction a large part of probability (around 
Wioc ~ 0.9) remains localized in the vicinity of the initial 
position of particles within the interval —5 < rc.1,2 < 5 
(Fig. 2). The numerical simulations show that in the in- 
teracting case the distribution mainly consists from two 
parts. One of them represents localized states and is 
frozen near the initial position of particles, another con- 
tinues to diffuse as in noninteractive case and corresponds 
to the Gaussian tail of distribution evolving in a diffusive 
way. Fig. 2 represents a typical distribution shape at an 
instant moment of time. These numerical data clearly 
demonstrate the qualitative change induced by interac- 
tion: appearence of localized component. 

If initially the particles are located on different sites 
then the value of Wi oc decreases with growth of the ini- 
tial distance Ar between them but its value still remains 
quite large if initial distance is about few sites. For fixed 
Ar the value of Wi oc is not sensitive to the initial choice 
of 774,2 and /9i,2- The dependence on time of probability 
to stay at the origin Pq averaged over the time interval 
[0, t] in shown in Fig. 3. Only in the noninteractive case 
Pq goes to zero with time while for nonzero U its value 
approaches to some constant. It is interesting to note 
that asymptotically Pq is larger than zero not only for 
A = 1 but also for A < 1 when the noninteracting case 
has continuous spectrum with waves ballistically propa- 
gating along the lattice. With the interaction decreas- 
ing the value of Pq decreases also but not very sharply 
(Fig. 3), also it is not sensitive to the sign of U. In the lo- 
calized phase A > 1 the interaction gives only a decrease 
of spreading over the lattice sites similarly to the case 
with A > 1. For example, in the case of Fig.l but with 
A = 1.05 and U = 1 the probability Pq is approximately 
9 times larger than for A = 1.05 and U = 0. 

The numerical results discussed above demonstrate 
that the TIP behavior in a quasiperiodic potential is quite 
different from the case of random potential where 
the interaction produces mainly delocalizing effect. The 
main reasons for this difference are probably the follow- 
ing. The derealization in the Harper model (1) at A = 1 
appears as the result of quantum tunneling between the 
sites with closes energies E n = 2A cos(?m + (3) which are 
exponentially far from each other but are close in en- 
ergy. Apparently, the interaction destroys these tiny res- 
onance conditions that leads to appearance of localized 
states. These states are localized in all directions on 2d 
lattice (711,712) (see Fig. 2). Therefore, they do not corre- 
spond to a situation in which interaction creats a coupled 
state which can propagate along the lattice. According 



to the numerical data these localized states are centered 
on the plane (774, 712) mainly along the diagonal n\ = ri2 
(two particles are close to each other) and their structure 
is approximately the same and independent of the posi- 
tion along the diagonal. This means that these states 
form a pure-point component in the energy spectrum. 
The question how this component is placed in respect to 
the spectrum of noninteracting particles remains open. 
One possibility is that this pure-point spectrum is lo- 
cated completely outside of noninteractive band [—8,8]. 
Such a case in some sense would be similar to an impu- 
rity state in a usual ballistic continuous band. Another 
possibility is that the pure-point spectrum is also par- 
tially located inside the band [—8, 8] in the gaps which 
exist for noninteractive problem. The second possibility 
looks to be more probable and more interesting. One of 
the indications in this direction is that strong decrease 
of U (Fig. 3) does not lead to disappearance of localized 
component. It is naturally to assume that the structure 
of the pure-point spectrum remains approximately the 
same for A < 1 (al least for not very small values of A). 
According to numerical data the singular-continuous part 
of the spectrum at A = 1 is not completely destroyed by 
the interaction, so that the quasi-diffusive spreading still 
takes place (Fig. 1), but it would be quite desirable to 
have rigorous mathematical results about the structure 
of the spectrum in the presence of interaction. 

The results discussed above were obtained for the TIP 
in the one-dimensional Harper model (3). They indi- 
cate that the interaction induced localized states also 
should exist in the original problem (2) of TIP on the 
two-dimensional lattice with magnetic flux. Indeed, here 
the interaction again should give a destruction of tiny 
resonance conditions for tunneling. However, the direct 
detailed investigations of the model (2) are required to 
make definite conclusions, but numerical studies of the 
model (2) are much more complicated than for the model 
(3). Finally, let us note that in the classical noninterac- 
tive limit the dynamics of models (2), (3) is integrable. 
Therefore, it would be interesting to study the effect of 
TIP in the kicked Harper model M where the one-particle 
classical dynamics is chaotic. 

In conclusion, the numerical investigations of TIP in 
the Harper model (3) show that the attractive/repulsive 
interaction leads to appearance of localized states and 
pure-point spectrum. This happens in the case when 
noninteractive system (U = 0) has quasi-diffusive wave 
packet spreading with singular-continuous spectrum (A = 
1) or even for A < 1 when it has the ballistic wave packet 
propagation and continuous spectrum. Such effect of in- 
teraction in quasiperiodic systems is attributed to the 
interaction induced destruction of tiny resonance con- 
ditions which in noninteractive system allowed to tun- 
nel between quasiresonant states leading to infinite wave 
packet spreading and decay of the probability to stay at 
the origin to zero. In the localized phase interaction gives 
only a decrease of localization length. Therefore, the sit- 
uation in the quasiperiodic systems is quite different from 
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the case of random potential where interaction between 
two particles gives an increase of localization length. 
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FIG. 1. Dependence of second moments 

o± =< (m±n2) 2 > on time t for TIP in the Harper model (3) 
(<t+ is full curve, ct_ is dashed). The parameters are A = 1, 
h = tt(5 1/2 - 1), (3i,2 = 2 1/2 ; U = for upper curves, U = 1 
for lower curves. System size isJVxiV = 301 x 301 sites, 
initially particles are at the same ni,2 = 0. 

FIG. 2. Dependence of integrated probability distribu- 
tions P + (full curve) and P- (dashed) on n± = (m ±n 2 ) 2 /2 
for the case of Fig.l and t = 4000. Left curves are for U = 1, 
right are for U = (shifted for clearity) . 



FIG. 3. Dependence of probability to stay at the initial 
state Po on time t for: A = 1, U = 1; A = 5/6,(7 = 5/6; 
A = 1, U = 1/4; A = 1, U = (curves from up to down) and 
7i,/?i,2 are as in Fig.l. 
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